Differential Gene Expression Analysis Research Articles

Carboxypeptidase M modulates BMSCs osteogenesis-adipogenesis via the MAPK/ERK pathway: Anintegrated single-cell and bulk transcriptomic study.

The pathogenesis of osteoporosis (OP) is closely associated with the disrupted balance between osteogenesis and adipogenesis in bone marrow-derived mesenchymal stem cells (BMSCs). We analyzed published single-cell RNA sequencing (scRNA-seq) data to dissect the transcriptomic profiles of bone marrow-derived cells in OP, reviewing 56 377 cells across eight scRNA-seq datasets from femoral heads (osteoporosis or osteopenia n = 5, osteoarthritis n = 3). Seventeen genes, including carboxypeptidase M (CPM), were identified as key osteogenesis-adipogenesis regulators through comprehensive gene set enrichment, differential expression, regulon activity, and pseudotime analyses. Invitro, CPM knockdown reduced osteogenesis and promoted adipogenesis in BMSCs, while adenovirus-mediated CPM overexpression had the reverse effects. Invivo, intraosseous injection of CPM-overexpressing BMSCs mitigated bone loss in ovariectomized mice. Integrated scRNA-seq and bulk RNA sequencing analyses provided insight into the MAPK/ERK pathway's role in the CPM-mediated regulation of BMSC osteogenesis and adipogenesis; specifically, CPM overexpression enhanced MAPK/ERK signaling and osteogenesis. In contrast, the ERK1/2 inhibitor binimetinib negated the effects of CPM overexpression. Overall, our findings identify CPM as a pivotal regulator of BMSC differentiation, which provides new clues for the mechanistic study of OP.

Upregulation of HMOX1 associated with M2 macrophage infiltration and ferroptosis in proliferative diabetic retinopathy

The roles of immune cell infiltration and ferroptosis in the progression of proliferative diabetic retinopathy (PDR) remain unclear. To identify upregulated molecules associated with immune infiltration and ferroptosis in PDR, GSE60436 and GSE102485 datasets were downloaded from the Gene Expression Omnibus (GEO). Genes associated with immune cell infiltration were examined through Weighted Gene Co-expression Network Analysis (WGCNA) and CIBERSORT algorithm. Common differentially expressed genes (DEGs) were intersected with ferroptosis-associated and immune cell infiltration-related genes. Localization of cellular expression was confirmed by single-cell analysis of GSE165784 dataset. Findings were validated by qRT-PCR, ELISA, Western blotting, and immunofluorescence staining. As a result, the infiltration of M2 macrophages was significantly elevated in fibrovascular membrane samples from PDR patients than the retinas of control subjects. Analysis of DEGs, M2 macrophage-related genes and ferroptosis-related genes identified three hub intersecting genes, TP53, HMOX1 and PPARA. qRT-PCR showed that HMOX1 was significantly higher in the oxygen-induced retinopathy (OIR) mouse model retinas than in controls. Single-cell analysis confirmed that HMOX1 was located in M2 macrophages. ELISA and western blotting revealed elevated levels of HMOX1 in the vitreous humor of PDR patients and OIR retinas, and immunofluorescence staining showed that HMOX1 co-localized with M2 macrophages in the retinas of OIR mice. This study offers novel insights into the mechanisms associated with immune cell infiltration and ferroptosis in PDR. HMOX1 expression correlated with M2 macrophage infiltration and ferroptosis, which may play a crucial role in PDR pathogenesis.

In vivo effect of vaginal seminal plasma application on the human endometrial transcriptome: a randomized controlled trial.

Studies in humans and animals suggest that seminal plasma, the acellular seminal fluid component, stimulates the endometrium to promote immune tolerance and facilitate implantation. We designed a randomized, double-blinded, placebo-controlled trial to investigate changes in the endometrial transcriptomic profile after vaginal application of seminal plasma. The study participants were randomized into two groups. Five women received a vaginal application of seminal plasma, and four received a placebo application with saline solution. The application was performed two days after human chorionic gonadotrophin-triggered ovulation in an unstimulated cycle. After 5-8 days, an endometrial biopsy was collected to analyze differences in the endometrial transcriptomic profile using microarray analyses. A differential gene expression analysis and a gene set analysis were performed. The gene set enrichment analysis showed a positive enrichment of pathways associated with the immune response, cell viability, proliferation and cellular movement. Moreover, pathways involved in implantation, embryo development, oocyte maturation, and angiogenesis were positively enriched. The differential gene expression analysis, after adjusting for multiple testing, showed no significantly differentially expressed genes between the two groups. A comparative analysis was also performed with similar studies conducted in other animals or in vitro using human endometrial cells. The comparative analysis showed that the effect of seminal plasma effect on the endometrium is similar in pigs, mice and in vitro human endometrial cells. The present study provides evidence that seminal plasma might impact the endometrium during the implantation window, with potential to affect endometrial receptivity and embryo development.

Large-scale transcriptomic and genomic analyses reveal a novel functional gene SERPINB6 for chicken carcass traits

BackgroundCarcass traits are crucial indicators of meat production efficiency. However, the molecular regulatory mechanisms associated with these traits remain unclear.ResultsIn this study, we conducted comprehensive transcriptomic and genomic analyses on 399 Tiannong partridge chickens to identify key genes and variants associated with carcass traits and to elucidate the underlying regulatory mechanisms. Based on association analyses with the elastic net (EN) model, we identified 12 candidate genes (AMY1A, AP3B2, CEBPG, EEF2, EIF4EBP1, FGFR1, FOXD3, GOLM1, LOC107052698, PABPC1, SERPINB6 and TBC1D16) for 4 carcass-related traits, namely live weight, dressed weight, eviscerated weight, and breast muscle weight. SERPINB6 was identified as the only overlapping gene by 3 analyses, EN model analysis, weighted gene co-expression network analysis and differential expression analysis. Cell-level experiments confirmed that SERPINB6 promotes the proliferation of chicken DF1 cells and primary myoblasts. Further expression genome-wide association study and association analysis indicated that rs317934171 is the critical site that enhances SERPINB6 expression. Furthermore, a dual-luciferase reporter assay proved that gga-miR-1615 targets the 3′UTR of SERPINB6.ConclusionsCollectively, our findings reveal that SERPINB6 serves as a novel gene for chicken carcass traits by promoting fibroblast and myoblast proliferation. Additionally, the downstream variant rs317934171 regulates SERPINB6 expression. These results identify a new target gene and molecular marker for the molecular mechanisms of chicken carcass traits.

Thymic gene expression analysis reveals a potential link between HIF-1A and Th17/Treg imbalance in thymoma associated myasthenia gravis

Myasthenia gravis (MG) is an immune-mediated disease frequently associated with thymic changes. Increased T helper 17 (Th17) cell activity and dysfunctional regulatory T (Treg) cells have been demonstrated in subgroups of MG. On the other hand, hypoxia-inducible factor 1 (HIF-1) has been shown to regulate the Th17/Treg balance by inducing Th17 differentiation while attenuating Treg development. To identify the underlying mechanisms of different thymic pathologies in MG development, we evaluated thymic samples from thymoma-associated myasthenia gravis (TAMG), MG with hyperplasia (TFH-MG) and thymoma without MG (TOMA) patients. Differential gene expression analysis revealed that TAMG and TFH-MG cells are associated with different functional pathways. A higher RORC/FOXP3 ratio provided evidence for Th17/Treg imbalance in TAMG potentially related to increased HIF1A. The hypoxic microenvironment in thymoma may be a driver of TAMG by increasing HIF1A. These findings may lead to new therapeutic approaches targeting HIF1A in the development of TAMG.

Characterization of alternative splicing events and prognostic signatures in gastric cancer

BackgroundAccumulating evidences indicate that the specific alternative splicing (AS) events are linked to the occurrence and prognosis of gastric cancer (GC). Nevertheless, the impact of AS is still unclear and needed to further elucidation.MethodsThe expression profile of GC and normal samples were downloaded from TCGA. AS events were achieved from SpliceSeq database. Cox regression together with LASSO analysis were employed to identify survival-associated AS events (SASEs) and calculate risk scores. PPI and pathway enrichment analysis were implemented to determine the function and pathways of these genes. Kaplan-Meier (K-M) analysis and Receiver Operating Characteristic Curves were used to evaluate the clinical significance of genes of SASEs. Q-PCR were applied to validate the hub genes on the survival prognosis in 47 GC samples. Drug sensitivity and immune cell infiltration analysis were conducted.ResultsIn total, 48 140 AS events in 10 610 genes from 361 GC and 31 normal samples were analyzed. Through univariate Cox regression, 855 SASEs in 763 genes were screened out. Further, these SASEs were analyzed by PPI and 17 hub genes were identified. Meanwhile, using Lasso and multivariate Cox regression analysis, 135 SASEs in 132 genes related to 7 AS forms were further screened and a GC prognostic model was constructed. K-M curves indicates that high-risk group has poorer prognosis. And the nomogram analysis on the basis of the multivariate Cox analysis was disclosed the interrelationships between 7 AS forms and clinical parameters in the model. Five key genes were then screened out by PPI analysis and Differential Expression Gene analysis based on TCGA and Combined-dataset, namely STAT3, RAD51B, SOCS2, POLE2 and TSR1. The expression levels of AS in STAT3, RAD51B, SOCS2, POLE2 and TSR1 were all significantly correlated with survival by qPCR verification. Nineteen drugs were sensitized to high-risk patients and eight immune cells showed significantly different infiltration between the STAD and normal groups.ConclusionsIn this research, the prognostic model constructed by SASEs can be applied to predict the prognosis of GC patients and the selected key genes are expected to become new biomarkers and therapeutical targets for GC treatment.

Particulate matter-induced gene expression patterns in human-derived cells based on 11 public gene expression datasets.

Exposure to particulate matter (PM) and house dust mite (HDM) can change the expression patterns of inflammation-, oxidative stress-, and cell death-related genes. We investigated the changes in gene expression patterns owing to PM exposure. This study examined the changes in gene expression patterns following PM exposure. We searched for differentially expressed genes (DEGs) following PM exposure using five cell line-based RNA-seq or microarray datasets and six human-derived datasets. The enrichment terms of the DEGs were assessed. DEG analysis yielded two gene sets. Thus, enrichment analysis was performed for each gene set, and the enrichment terms related to respiratory diseases were presented. The intersection of six human-derived datasets and two gene sets was obtained, and the expression patterns following PM exposure were observed. Two gene sets were obtained for cells treated with PM and their expression patterns were presented following verification in human-derived cells. Our findings suggest that exposure to PM2.5 and HDM may reveal changes in genes that are associated with diseases, such as allergies, highlighting the importance of mitigating PM2.5 and HDM exposure for disease prevention.

Identification of driving genes of familial adenomatous polyposis by differential gene expression analysis and weighted gene co-expression network analysis.

Despite the advancement of new screening strategies and the advances in pharmacological therapies, the cancerization rates of familial adenomatous polyposis (FAP) are stable and even increased in the last years. Therefore, it necessitates additional research to characterize and understand the underlying mechanisms of FAP. To determine the genes that drive the pathogenesis of familial adenomatous polyposis (FAP). We performed on a cohort (GSE111156) gene profile, which consist of four group of gene expressions (the gene expressions of cancer, adenoma and normal tissue of duodenal cancer from patients with FAP were defined as Case N, Case A and Case C respectively, while that of adenoma tissue from patients with FAP who did not have duodenal cancer was Ctrl A). Tracking Tumor Immunophenotype (TIP) website was applied to reveal immune infiltration profile and signature genes of FAP. We merged the genes of key module (pink and midnight module) with signature genes to obtained the biomarkers related with FAP pathogenesis. The expression of these five biomarkers in FAP intratumoral region (IT) and tumor rim (TR) was detected with Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR). In total, 220, 23 and 63 DEGs were determined in Cases C, A and N, in comparison to Ctrl A. In total, 196 and 10 DEGs were determined in Cases C and A, separately, as compared to Case N. A total of four biomarkers including CCL5, CD3G, CD2 and TLR3 were finally identified associated with pink module, while only one biomarker (KLF2) associated with midnight module was identified. All biomarkers were evidently raised in FAP IT tissues utilizing qRT-PCR. We identified five potential biomarkers for pathogenesis of FAP to understand the fundamental mechanisms of FAP progression and revealed some probable targets for the diagnosis or treatment of FAP.

Integrated time-series biochemical, transcriptomic, and metabolomic analyses reveal key metabolites and signaling pathways in the liver of the Chinese soft-shelled turtle (Pelodiscus sinensis) against Aeromonas hydrophila infection

IntroductionAeromonas hydrophila, a bacterium widely distributed in the natural environment, causes multiple diseases in various animals. Exploring the mechanism of the host defense against A. hydrophila can help develop efficient strategies against Aeromonas infection.MethodsHerein, we investigated the temporal influence of A. hydrophila on the Chinese soft-shelled turtle, an economically important species, at the biochemical, transcriptomic, and metabolomic levels. Plasma parameters were detected with the test kits. Transcriptome and metabolome were respectively applied to screen the differentially expressed genes and metabolites.ResultsThe contents or activities of these plasma parameters were significantly increased at 24 hpi and declined at 96 hpi, indicating that 24 and 96 hpi were two important time points during infection. Totals of 3121 and 274 differentially expressed genes (DEGs) from the transcriptome while 74 and 91 differentially abundant metabolites (DAMs) from the metabolome were detected at 24 and 96 hpi. The top DEGs at 24 hpi included Ccl2, Ccl3, Ccl4, Il1β, Il6, Il7, Il15, Tnf, and Tnfr1 while Zap70, Cd3g, Cd8a, Itk, Pik3r3, Cd247, Malt1, and Cd4 were the most abundant at 96 hpi. The predominant DAMs included O-phospho-L-serine, γ-Aminobutyric acid, orotate, L-tyrosine, and L-tryptophan at 24 hpi, as well as L-glutamic acid, L-arginine, glutathione, glutathione disulfide, and citric acid at 96 hpi.DiscussionThe combined analysis of DEGs and DAMs revealed that tryptophan metabolism, nicotinate and nicotinamide metabolism, as well as starch and sucrose metabolism, were the most important signaling pathways at the early infective stage while tyrosine metabolism, pyrimidine metabolism, as well as alanine, aspartate and glutamate metabolism were the most crucial pathways at the later stage. In general, our results indicated that the Chinese soft-shelled turtle displays stage-specific physiological responses to resist A. hydrophila infection.

Identification of key genes affecting ventilator-induced diaphragmatic dysfunction in diabetic mice

BackgroundMechanical ventilation (MV) is often required in critically ill patients. However, prolonged mechanical ventilation can lead to Ventilator-induced diaphragmatic dysfunction (VIDD), resulting in difficulty in extubation after tracheal intubation, prolonged ICU stay, and increased mortality. At present, the incidence of diabetes is high in the world, and the prognosis of diabetic patients with mechanical ventilation is generally poor. Therefore, the role of diabetes in the development of VIDD needs to be discovered.MethodsMV modeling was performed on C57 mice and DB mice, and the control group was set up in each group. After 12 h of mechanical ventilation, the muscle strength of the diaphragm was measured, and the muscle fiber immunofluorescence staining was used to verify the successful establishment of the MV model. RNA sequencing (RNA-seq) method was used to detect mRNA expression levels of the diaphragms of each group, and then differential expressed gene analysis, Heatmap analysis, WGCNA analysis, Venn analysis, GO and KEGG enrichment analysis were performed. qRT-PCR was used to verify the expression of the selected mRNAs.ResultsOur results showed that, compared with C57 control mice, the muscle strength and muscle fiber cross-sectional area of mice after mechanical ventilation decreased, and DB mice showed more obvious in this respect. RNA-seq showed that these differential expressed (DE) mRNAs were mainly related to genes such as extracellular matrix, collagen, elastic fiber and Fbxo32. GO and KEGG enrichment analysis showed that the signaling pathways associated with diabetes were mainly as follows: extracellular matrix (ECM), protein digestion and absorption, PI3K-Akt signaling pathway, calcium signaling pathway, MAPK signaling pathway and AGE-RAGE signaling pathway in diabetic complications, etc. ECM has the closest relationship with VIDD in diabetic mice. The key genes determined by WGCNA and Venn analysis were validated by quantitative real-time polymerase chain reaction (qRT-PCR), which exhibited trends similar to those observed by RNA-seq.ConclusionVIDD can be aggravated in diabetic environment. This study provides new evidence for mRNA changes after mechanical ventilation in diabetic mice, suggesting that ECM and collagen may play an important role in the pathophysiological mechanism and progression of VIDD in diabetic mice, and provides some clues for the research, diagnosis, and treatment of VIDD in diabetic context.

Functional transcriptome analysis revealed upregulation of MAPK-SMAD signalling pathways in chronic heat stress in crossbred cattle.

Animal geneticists and breeders have the impending challenge of enhancing the resilience of Indian livestock to heat stress through better selection strategies. Climate change's impact on livestock is more intense in tropical countries like India where dairy cattle crossbreeds are more sensitive to heat stress. The main reason for this study was to find the missing relative changes in transcript levels in thermo-neutral and heat stress conditions in crossbred cattle through whole-transcriptome analysis of RNA-Seq data. Differentially expressed genes (DEGs) identified based on the minimum log twofold change value and false discovery rate 0.05 revealed 468 up-regulated genes and 2273 down-regulated significant genes. Functional annotation and pathway analysis of these significant DEGs were compared based on Gene Ontology (Biological process), Kyoto Encyclopedia of Genes and Genome (KEGG), and Reactome pathways using g: Profiler, ShinyGO v0.76, and iDEP.951 web tools. On finding network visualization, the most over-represented and correlated pathways were neuronal and sensory organ development, calcium signalling pathway, Mitogen-activated protein kinase (MAPK) and Smad signalling pathway, Ras-proximate-1, or Ras-related protein 1 (Rap 1) signalling pathway, apoptosis, and oxidative stress. Similarly, down-regulated genes were most expressed in mRNA processing, immune system, B-cell receptor signalling pathway, Nucleotide oligomerization domain(NOD)-like receptors (NLRs) signalling pathway and nonsense-mediated decay (NMD) pathway. The heat stress-responsive genes identified in this study will facilitate our understanding of the molecular basis for climate resilience and heat tolerance in Indian dairy crossbreeds.

Transcriptomics reveals age-related changes in ion transport–related factors in yak lungs

Yaks inhabit high-altitude, low-oxygen regions, where ion transport functions play a crucial role in maintaining intracellular and extracellular ionic balance and regulating pulmonary vascular tension. These functions affect pulmonary ventilation and blood flow rate, aiding tissue development and enhancing oxygen transfer efficiency, thus facilitating better adaptation to hypoxic environments. To investigate the regulatory mechanisms of ion transport-related factors on the growth and development of yak lungs, we employed RNA sequencing (RNA-seq)for sequencing the transcriptome in the lung tissues of neonatal (1-day-old), juvenile (1-year-old), and adult (4-year-old) yaks. We also performed differential gene expression and functional analyses. The results yielded 26 genes associated with ion transport, mainly enriched in the salivary and pancreatic secretion pathways. Finally, we used several methods including quantitative polymerase chain reaction (qRT-PCR), and Western blotting (WB), immunohistochemical (IHC) and immunofluorescence (IF) staining to determine the distribution of the expression of the ion transport genes FOXI1, KCNMA1, and SLC12A2 in yak lung tissues. qRT-PCR and WB results indicated that mRNA and protein relative expression levels of FOXI1 and SLC12A2 were significantly higher in neonatal yaks than in juvenile and adult yaks (all p < 0.05), whereas those of KCNMA1 were significantly higher in adult yaks than in neonatal and juvenile yaks (all p < 0.05). IHC and IF results demonstrated that FOXI1, KCNMA1, and SLC12A2 were distributed among the epithelial mucosal layers (including ciliated, goblet, and Clara cells) of the yaks’ bronchi and their branches in the lungs across different age groups of yak. Therefore, our results suggested that FOXI1, KCNMA1, and SLC12A2 may be strongly associated with the development and aging processes in yak lungs. These results provide insights into the molecular mechanisms underlying the yak’s adaptation to high-altitude environments and valuable references for further research.

Matrisome Transcriptome Dynamics during Tissue Aging

The extracellular matrix (ECM) is a complex three-dimensional network of macromolecules that provides structural support for the cells and plays a significant role in tissue homeostasis and repair. Growing evidence indicates that dysregulation of ECM remodeling contributes to various pathological conditions in the body, including age-associated diseases. In this work, gene expression data of normal human tissues obtained from the Genotype-Tissue Expression project, as well as data from MatrisomeDB 2.0, the ECM-protein knowledge database, are used to estimate the age-dependent matrisome transcriptome dynamics in the blood, heart, brain, liver, kidneys, lungs, and muscle. Differential gene expression (DE) analysis revealed dozens of matrisome genes encoding both structural elements of the ECM and ECM-associated proteins, which had a tissue-specific expression profile with age. Among common DE genes that changed expression with age in at least three tissues, COL18A1, MFAP1, IGFBP7, AEBP1, LTBP2, LTBP4, LG14, EFEMP1, PRELP, BGN, FAM20B, CTSC, CTSS, and CLEC2B were observed. The findings of the study also reveal that there are sex-specific alterations during aging in the matrisome gene expression. Taken together, the results obtained in this work may help in understanding the role of the ECM in tissue aging and might prove valuable for the future development of the field of ECM research in general.

Identification of Prognostic Genes in Acute Myeloid Leukemia Microenvironment: A Bioinformatic and Experimental Analysis.

Acute myeloid leukemia (AML) is a lethal hematologic malignancy with a variable prognosis that is highly dependent on the bone marrow microenvironment. Consequently, a better understanding of the AML microenvironment is crucial for early diagnosis, risk stratification, and personalized therapy. In recent years, the role of bioinformatics as a powerful tool in clarifying the complexities of cancer has become more prominent. Gene expression profile and clinical data of 173 AML patients were downloaded from the TCGA database, and the xCell algorithm was applied to calculate the microenvironment score (MS). Then, the correlation of MS with FAB classification, and CALGB cytogenetic risk category was investigated. Differentially expressed genes (DEGs) were identified, and the correlation analysis of DEGs with patient survival was done using univariate cox. The prognostic value of candidate prognostic DEGs was confirmed based on the GEO database. In the last step, real-time PCR was used to compare the expression of the top three prognostic genes between patients and the control group. During TCGA data analysis, 716 DEGs were identified, and survival analysis results showed that 152 DEGs had survival-related changes. In addition, the prognostic value of 31 candidate prognostic genes was confirmed by GEO data analysis. Finally, the expression analysis of FLVCR2, SMO, and CREB5 genes, the most related genes to patients' survival, was significantly different between patients and control groups. In summary, we identified key microenvironment-related genes that influence the survival of AML patients and may serve as prognostic and therapeutic targets.

Deciphering novel TCF4-driven mechanisms underlying a common triplet repeat expansion-mediated disease.

Fuchs endothelial corneal dystrophy (FECD) is an age-related cause of vision loss, and the most common repeat expansion-mediated disease in humans characterised to date. Up to 80% of European FECD cases have been attributed to expansion of a non-coding CTG repeat element (termed CTG18.1) located within the ubiquitously expressed transcription factor encoding gene, TCF4. The non-coding nature of the repeat and the transcriptomic complexity of TCF4 have made it extremely challenging to experimentally decipher the molecular mechanisms underlying this disease. Here we comprehensively describe CTG18.1 expansion-driven molecular components of disease within primary patient-derived corneal endothelial cells (CECs), generated from a large cohort of individuals with CTG18.1-expanded (Exp+) and CTG 18.1-independent (Exp-) FECD. We employ long-read, short-read, and spatial transcriptomic techniques to interrogate expansion-specific transcriptomic biomarkers. Interrogation of long-read sequencing and alternative splicing analysis of short-read transcriptomic data together reveals the global extent of altered splicing occurring within Exp+ FECD, and unique transcripts associated with CTG18.1-expansions. Similarly, differential gene expression analysis highlights the total transcriptomic consequences of Exp+ FECD within CECs. Furthermore, differential exon usage, pathway enrichment and spatial transcriptomics reveal TCF4 isoform ratio skewing solely in Exp+ FECD with potential downstream functional consequences. Lastly, exome data from 134 Exp- FECD cases identified rare (minor allele frequency <0.005) and potentially deleterious (CADD>15) TCF4 variants in 7/134 FECD Exp- cases, suggesting that TCF4 variants independent of CTG18.1 may increase FECD risk. In summary, our study supports the hypothesis that at least two distinct pathogenic mechanisms, RNA toxicity and TCF4 isoform-specific dysregulation, both underpin the pathophysiology of FECD. We anticipate these data will inform and guide the development of translational interventions for this common triplet-repeat mediated disease.

Transcriptomic responses to drought stress in the Patagonian southern beech Nothofagus alpina

BackgroundDeciphering the genetic architecture of drought tolerance could allow the candidate genes identification responding to water stress. In the Andean Patagonian forest, the genus Nothofagus represents an ecologically relevant species to be included in different genomic studies. These studies are scarce in South American ecosystems however represent an important source of genomic data in order to interpret future climate-change environment scenarios of these emblematic forests. Here, we achieved the assemblage of the transcriptome of N. alpina while searching for key genes of activated or suppressed metabolic pathways in response to drought stress.ResultsDe novo transcriptome assembly resulted in 104,030 transcripts. Following confirmation of drought conditions, based on reduction of leaf water potential and stomatal conductance, a differential gene expression analysis resulted in 2720 significantly expressed genes (1601 up-regulated and 1119 down-regulated). Enrichment analysis (over-representation analysis and gene set enrichment analysis) resulted in more than one hundred stress-responsive term ontologies (i.e. biological processes) and pathways. Terms such as response to abscisic acid and pathways such as plant hormone signal transduction or starch and sucrose metabolism were over-represented. Protein–protein interaction assessment resulted in networks with significantly expressed top common hub gene clusters (e.g. plant-type cell wall biogenesis among down-regulated or ABA-signalling among up-regulated). These networks evidenced important regulators at gene expression such as transcriptional factors.ConclusionsResponses of N. alpina seedlings to drought stress were evidenced by the activation of several genes linked to GO biological processes and KEGG pathways, which were mainly based on over-expression of specific protein kinases, phosphatases, synthases and transcription factors. This suggests an up-regulation of signalling pathways, triggered through plant hormones such as abscisic acid or auxin, which could counteract the osmotic stress created as a probable immediate response to drought. On the other hand, groups of carbon fixation genes related to the galactose metabolism, photosynthesis, secondary wall biogenesis, and fatty acid biosynthesis degradation were down-regulated under drought. Overall, our results provide new genomic data for understanding how non-model long-lived trees of Patagonian forests would acclimate to environmental changes.

Heterozygous RB1 mutation enhanced ATP production in human iPSC-derived retinal organoids.

Recent in vitro studies using RB1+/- fibroblasts and MSCs have shown molecular and functional disruptions without the need for biallelic loss of RB1. However, this was not reflected in the recent in vitro studies employing RB1+/- retinal organoids. To gain further insights into the molecular disruptions in the RB1+/- retinal organoids, we performed a high throughput RNA sequencing analysis. iPSCs were generated from RB1+/+ and RB1+/- OAMSCs derived from retinoblastoma patients. RB1+/+ and RB1+/- iPSCs were subjected to a step-wise retinal differentiation protocol. Retinal differentiation was evaluated by Real-time PCR and flow cytometry analysis of the retinal markers. To gain further insights into the molecular differences in RB1+/- retinal organoids, a high throughput RNA sequencing followed by differential gene expression analysis and gene set enrichment analysis (GSEA) was performed. The analysis revealed a shift from the regular metabolic process of glycolysis to oxidative phosphorylation in the RB1+/- retinal organoids. To investigate further, we performed assays to determine the levels of pyruvate, lactate and ATP in the retinal organoids. The results revealed significant increase in ATP and pyruvate levels in RB1+/- retinal organoids of day 120 compared to that of the RB1+/+. The results thus revealed enhanced ATP production in the RB1+/- retinal organoids. The study provides novel insights into the metabolic phenotype of heterozygous RB1 mutant suggesting dysregulation of energy metabolism and glycolytic pathways to be first step even before the changes in cellular proliferation or other phenotypic consequences ensue.

Transcriptome profiles of Trypanosoma brucei rhodesiense in Malawi reveal focus specific gene expression profiles associated with pathology.

Sleeping sickness caused by Trypanosoma brucei rhodesiense is a fatal disease and endemic in Southern and Eastern Africa. There is an urgent need to develop novel diagnostic and control tools to achieve elimination of rhodesiense sleeping sickness which might be achieved through a better understanding of trypanosome gene expression and genetics using endemic isolates. Here, we describe transcriptome profiles and population structure of endemic T. b. rhodesiense isolates in human blood in Malawi. Blood samples of r-HAT cases from Nkhotakota and Rumphi foci were collected in PaxGene tubes for RNA extraction before initiation of r-HAT treatment. 100 million reads were obtained per sample, reads were initially mapped to the human genome reference GRCh38 using HiSat2 and then the unmapped reads were mapped against Trypanosoma brucei reference transcriptome (TriTrypDB54_TbruceiTREU927) using HiSat2. Differential gene expression analysis was done using the DeSeq2 package in R. SNP calling from reads that were mapped to the T. brucei genome was done using GATK in order to identify T.b. rhodesiense population structure. 24 samples were collected from r-HAT cases of which 8 were from Rumphi and 16 from Nkhotakota foci. The isolates from Nkhotakota were enriched with transcripts for cell cycle arrest and stumpy form markers, whereas isolates in Rumphi focus were enriched with transcripts for folate biosynthesis and antigenic variation pathways. These parasite focus-specific transcriptome profiles are consistent with the more virulent disease observed in Rumphi and a less symptomatic disease in Nkhotakota associated with the non-dividing stumpy form. Interestingly, the Malawi T.b. rhodesiense isolates expressed genes enriched for reduced cell proliferation compared to the Uganda T.b. rhodesiense isolates. PCA analysis using SNPs called from the RNAseq data showed that T. b. rhodesiense parasites from Nkhotakota are genetically distinct from those collected in Rumphi. Our results suggest that the differences in disease presentation in the two foci is mainly driven by genetic differences in the parasites in the two major endemic foci of Rumphi and Nkhotakota rather than differences in the environment or host response.

Integrated anatomical structure, physiological, and transcriptomic analyses to identify differential cold tolerance responses of Ziziphus jujuba mill. ‘Yueguang’ and its autotetraploid ‘Hongguang’

Cold stress is a limiting stress factor that limits plant distribution and development; however, polyploid plants have specific characteristics such as higher resistance to abiotic stress, especially cold stress, that allow them to overcome this challenge. The cultivated cultivar Ziziphus jujuba Mill. ‘Yueguang’ (YG) and its autotetraploid counterpart ‘Hongguang’ (HG) exhibit differential cold tolerance. However, the underlying molecular mechanism and methods to enhance their cold tolerance remain unknown. Anatomical structure and physiological analysis indicated YG had a higher wood bark ratio, and xylem ratio under cold treatment compared to HG. However, the half-lethal temperature (LT50), cortex ratio, and malondialdehyde (MDA) content were significantly decreased in YG than HG, which indicated YG was cold tolerant than HG. Transcriptome analysis showed that 2084, 1725, 2888, and 2934 differentially expressed genes (DEGs) were identified in HC vs YC, H20 vs Y20, Y20 vs YC, and H20 vs HC treatment, respectively. Meanwhile, KEGG enrichment analysis of DEGs showed that several metabolic pathways, primarily plant hormone signal transduction and the MAPK signaling pathway, were involved in the differential regulation of cold tolerance between YG and HG. Furthermore, exogenous abscisic acid (ABA) and brassinolide (BR) treatments could improve their cold tolerance through increased SOD and POD activities, decreased relative electrical conductivity, and MDA content. All of these findings suggested that plant hormone signal transduction, particularly ABA and BR, might have an important role in the regulation of differential cold tolerance between YG and HG, laying the foundation for further improving cold tolerance in jujube and examining the molecular mechanisms underlying differences in cold tolerance among different ploidy cultivars.

Deciphering Müller cell heterogeneity signatures in diabetic retinopathy across species: an integrative single-cell analysis

Diabetic retinopathy (DR), a leading cause of visual impairment, demands a profound comprehension of its cellular mechanisms to formulate effective therapeutic strategies. Our study presentes a comprehensive single-cell analysis elucidating the intricate landscape of Müller cells within DR, emphasizing their nuanced involvement. Utilizing scRNA-seq data from both Sprague–Dawley rat models and human patients, we delineated distinct Müller cell clusters and their corresponding gene expression profiles. These findings were further validated through differential gene expression analysis utilizing human transcriptomic data. Notably, certain Müller cell clusters displayed upregulation of the Rho gene, implying a phagocytic response to damaged photoreceptors within the DR microenvironment. This phenomenon was consistently observed across species. Additionally, the co-expression patterns of RHO and PDE6G within Müller cell clusters provided compelling evidence supporting their potential role in maintaining retinal integrity during DR. Our results offer novel insights into the cellular dynamics of DR and underscore Müller cells as promising therapeutic targets for preserving vision in retinal disorders induced by diabetes.